Lagrangian statistics and coherent structures in two-dimensional turbulence 
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Measurements of Lagrangian single-point and multiple-point statistics in a quasi-two-dimensional 
stratifed layer system are reported. The system consists of a layer of salt water over an immiscible 
layer of Fluorinert and is forced electromagnetically so that mean-squared vorticity is injected at a 
well-defined spatial scale. Simultaneous cascades develop in which enstrophy flows predominately to 
small scales whereas energy cascades, on average, to larger scales. Lagrangian velocity correlations 
and one- and two-point displacements are measured for random initial conditions and for initial 
positions within topological centers and saddles. The behavior of these quantities can be understood 
in terms of the trapping characteristics of long-lived centers, the slow motion near strong saddles, 
and the rapid fluctuations outside of either centers or saddles. 



I. INTRODUCTION 

There is a growing amount of experimental and numer- 
ical evidence indicating that two-dimensional (2D) turbu- 
lence at moderate to low Reynolds number is dominated 
by the existence of coherent structures|l|,|2j. This obser- 
vation is more or less independent of whether the turbu- 
lence maintains an energetically steady state with a fully 
developed inverse energy cascade or merely coarsens with 
time as in decaying 2D turbulence. Coherent structure 
dominance and the unique aspects of the 2D dual cas- 
cade picture, i.e. an inverse cascade moving energy from 
small to large length scales and a direct cascade moving 
mean-squared vorticity, or enstrophy, from large to small 
length scales, clearly distinguishes 2D turbulence from 
its three-dimensional (3D) counterpart. An emphasis in 
the study of 2D turbulence is to understand what im- 
pact these structures have on the statistics of turbulence. 
Recent numerical and experimental results demonstrate 
that statistical quantities display different behavior in 
regions around locally hyperbolic points in the velocity 
field (saddles) and within vortices (centers) 0, 3 ■ 

The existence of coherent structures is expected to play 
an important role in the turbulent dispersion (mixing) of 
passive scalarsjj^. In particular, centers, if long lived, 
should act as trapping regions that prevent scalar advec- 
tion over their encircling manifolds. On the other hand, 
saddles should quite effectively enhance the dispersion by 
the stretching of fluid elements into thin filaments along 
the saddles unstable manifold. It is this picture that un- 
derlies the description of trajectories of fluid elements in 
2D turbulence in terms of successive trapping and flight 
events, i.e., Levy flightsjfj. Although there is little ev- 
idence supporting this idea, the proposal highlights the 
importance of accounting for coherent structures in any 
mixing model. 

How coherent structures affect the mixing properties 
of 2D turbulent flows can have profound impact on 
real world problems. An example of geophysical signif- 
icance is the Antarctic circumpolar vortex which, dur- 
ing the winter months, prevents ozone-depleted air over 
the South Pole from mixing with ozone-rich air at mid 



latitudes^- Also, the possibility that 3D turbulence may 
not be as dominated by coherent structures as is 2D tur- 
bulence does not mean that coherent structures do not 
play an important role in 3D mixingj^]. Understanding 
how coherent structures influence mixing in 2D may aid 
in the understanding of 3D mixing. 

The mixing problem is generally approached by ap- 
plying statistical treatments to fluid-element (or parti- 
cle) trajectories, i.e., to the Lagrangian dynamics of the 
fluid. Classical descriptions of Lagrangian dynamics in- 
volve the statistical properties of both single trajectories 
(N = 1) following the work of G.I. Taylor|9j and relative 
displacement of two trajectories (N — 2) following the 
work of Richardson [Toj. In this manuscript, the effects of 
coherent structures on mixing in 2D turbulence will be 
quantified by measuring their effects on these classical 
quantities. Recently there has been important theoreti- 
cal and numerical work on the evolution of fluid patches 
(groups of N fluid elements with N > 2) 11]. This ap- 
proach is interesting and could include the effects of co- 
herent structures on the evolution of fluid patches but is 
beyond the scope of the present work. 

Following the discussion in the classical descrip- 
tions of particle motion in 2D fluids start from the obser- 
vation that the Lagrangian trajectory, r = r(t), of a fluid 
element or non-inertial particle embedded in the fluid is 
governed by the stochastic equation: 



dr = u(r ,t)dt + 2(2K)*/3(i). 



(1) 



In Eq. u(r, t) is the incompressible fluid velocity at 
position r at time t, K is the molecular diffusivity and f3 
represents standard Brownian motion. Since our system 
is in a regime dominated by turbulent advection and not 
by thermal fluctuations, k will be set to zero. 

For single-particle trajectories, the important statisti- 
cal property is the distance a particle has traveled, on 
average, over a given time interval St = t — to. Defining 
Sr = r(t) — r(io) one can show that 



±(Sr) 2 = 2u(r ■(*),*). /* u(r(s), s)ds, 



(2) 



which, upon performing an ensemble average and assum- 
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ing u is statistically stationary, becomes an equation for 
the second moment of the single particle displacement 



> )2 > = 2 /> 



(r(to),to).u(r(a),a))ds. (3) 



The behavior of ((Sr) 2 ) depends on how St compares with 
the Lagrangian correlation time t£ = J t (u(r(to)) to) ■ 
u(r(s),s))ds/(u 2 ). If St < t l then u(R(i),*) « 
u(r(to)j so that the term within the integral of Eq. [31 
simply becomes (u 2 ). In this case ((Sr) 2 ) sa (u 2 )St 2 , i.e. 
ballistic motion. If, on the other hand, St 3> tl then 
the integral in Eq. [3| reduces to tx(u 2 ) yielding Brow- 
nian walk behavior for the single particle displacement: 
((Sr) 2 ) w 2(u 2 )T£(5i. These limits follow simply from the 
properties of the velocity autocorrelation. For intermedi- 
ate times the behavior of ((<5r) 2 ) depends on the nature of 
the inertial turbulent fluctuations and not much is known 
about its scaling in 2D turbulence^]- 

For two-particle trajectories, the quantity of interest is 
the statistics of the relative displacement of particles in a 
turbulent flow. Defining the relative displacement of two 
trajectories R12 = ^(t) — ri(t), Eq. 2] implies that 



dt 



R12 = u(r 2 ,t) - u(r l5 t) = Au. 



(4) 



As in the case of the single trajectories, an equation for 
R 2 can be obtained (dropping the subscripts for simplic- 
ity): 

dR 2 



dt 



= 2R • Au. 



(5) 



If R is small, Au should be linearly proportional to R, 
i.e., the first term in a Taylor series. This linear behavior 
yields exponential growth in time for R 2 . On the other 
hand, if R is within the inertial energy range one expects 



2/3 



St 



that Au • R w IR 1 / 3 ^. This yields |R| 2 / 3 - R\ 
where Rq is the trajectory spacing at time to- Under the 
assumption that Rq is small compared to the scales of 
interest we get |R| 2 « (t — to) 3 , the famous t 3 law derived 
by Richardson[ll|. A Richardson range was reported in 
both 2D 01 an d 3D experimental systems (l^. 

Our work on 2D Lagrangian statistics and coherent 
structures begins with a description and characterization 
of the experimental apparatus and measurement tech- 
niques in Sec. [H] Section IIIII presents the basic single 
trajectory and two trajectory Lagrangian statistics. In 
Sec. IIVI the quantities presented in Sec. IIIII are condi- 
tioned on the presence of coherent structures to elucidate 
the role that such structures play in the mixing process. 
Finally, in Sec. [3 the implications of and future opportu- 
nities for these Lagrangian measurements of turbulence 
in 2D systems are discussed. 



II. EXPERIMENTAL METHODS 



of experiments, as well as the technique used for obtain- 
ing velocity and trajectory information, namely particle 
tracking. In Sec. Ill Bl the technique for identifying co- 
herent structures within the flow is established. Finally, 
Sec. Ill CI presents basic statistical properties of the tur- 
bulence produced in the layer. 



A. System and measurement 

The two-point prediction of (|R| 2 ) ~ t 3 obtained in 
Sec. Q] requires the presence of an inverse energy cascade 
where (Su 2 ^) scales as R 2 ! 3 (as we will see this is not a 
strict requirement). For 2D, or more pointedly quasi-2D, 
turbulence this implies that the turbulence must be con- 
tinuously forced. One method for continuously forcing 
the body of a fluid, originally pioneered by Dolzhansky in 
1979 |15| . is to subject a current carrying fluid to external 
magnetic fields. This method of forcing has since evolved 
into the stratified electromagnetic layer [l^ which has be- 
come one of the more common systems for the study of 
2D turbulence. 

The stratified electromagnetic layer is comprised of a 
dense salt water layer, typically of order 0.3 cm deep and 
20 cm x 20 cm square, underneath a less dense solution 
of roughly the same depth. A current is passed in the 
plane of the layer and the layers are subject to a spatially 
varying magnetic field penetrating them vertically. The 
resultant Lorentz force drives fluid motion. The evolution 
of the stratified layer system is expected to approximate 
the forced/damped 2D Navier-Stokes equation, 



duj 
dt 



' dx s 



dp 

dxi 



dx 2 



+ F Z . 



(6) 



supplemented by the incompressibility condition diUi = 
0. As usual, u is the fluid velocity field, p is the den- 
sity normalize pressure, F is the external electromagnetic 
forcing, v is the fluid kinematic viscosity, and Einstein 
summation is used throughout. The linear term with co- 
efficient a represents the effects of frictional drag owing 
to the container bottom. 

In the set of experiments presented here, the system de- 
scribed above is modified slightly by replacing the lower 
layer of fluid with Fluorinert FC-75 and the upper layer 



contacts 



+ 



salt water 

4= 



fluorinert ip^i 



magnet array 



Section lll Al describes the stratified electromagnetically 
forced layer that was the quasi 2D fluid used in this set 



FIG. 1: Schematic illustration of the stratified layer 2D tur- 
bulence experiment. See Sec. Ill Al for detailed description. 
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FIG. 2: Raw particle tracks obtained from the stratified layer 
over four consecutive frames (0.07 s). 

by a dense salt water solution of 20% by mass NaCl with 
a small amount of liquid detergent added to lower surface 
tension and help with dissolution of tracer particles. The 
Fluorinert is used because it has a density 1.8 times that 
of water with near the same viscosity, which allows for 
much stronger stratification than in the case of two salt 
water solutions. It is also a strong dielectric so the lower 
layer is completely passive and only the upper salt water 
layer is electromagnetically forced. Finally, and perhaps 
most importantly, the Fluorinert and water are immisci- 
ble. These combined features allow the Fluorinert system 
to maintain stratification indefinitely, allowing the salt 
water layer to be driven harder than in previous systems 

The experimental apparatus, shown schematically in 
Fig. ^ consisted of a 0.3 cm thick layer of salt-solution 
suspended over a 0.3 cm thick layer of Fluorinert. The 
layers were contained in a 40 cm x 20 cm box with reser- 
voirs at each end across which copper electrodes were 
placed in the fluid. For the results reported in this paper 
an alternating square- wave current with frequency 0.5 Hz 
and amplitude 0.75 amps was driven through the salt so- 
lution. A set of 1.27 cm diameter rare-earth magnets of 
approximately 0.7 T residual field strength were arranged 
with alternating field direction in a 20 cm x 20 cm square 
array with a period of 2.54 cm and oriented at 45° with 
respect to the current direction. The combination of the 
current and the magnetic field produces a Kolmogorov- 
like forcing 0, of alternating shear bands with the 
shear direction along y and periodicity r^j rj 1.8 cm 
in the x direction (which implies hi n j = 27r/ri n j = 3.5 
rad/cm). Using the layer depths of 0.3 cm yields an ap- 
proximation of a w 0.125 Hz for the frictional coupling 
assuming a simple linear shear in the Fluorinert. Also, 
the salt solution upper layer has a viscosity around 1.15 
that of water. 



To obtain velocity information as well as trajectory in- 
formation used in generating Lagrangian statistics, the 
upper salt-solution layer was seeded with polycrystalline 
powder with mean diameter 75 fim and density 0.98 
gm/cc. Images of the particle fields, illuminated using 
several Xenon short-arc flash lamps, were obtained with 
a 1280 x 1024 pixel CCD camera at a frame rate of 60 
Hz. The velocity field was obtained from image pairs us- 
ing particle tracking velocimetry derived from two earlier 
methods fl8l ITgj . For a typical pair of images, of order 
2.5 x 10 4 particle tracks were obtained which were inter- 
polated to a 126 x 100 velocity field array. A typical raw 
particle track field is shown in Fig. 

From this raw data one can obtain the required particle 
trajectories in two ways: splice together the raw particle 
tracks or use the interpolated Eulerian fields and generate 
particle trajectories by solving the advection equation. 
Both techniques will be used here. Most of the results 
will rely, however, on the generated tracks because it is 
difficult to obtain enough statistics with particle track 
splicing (an error rate of only 1 percent splicing from one 
frame to the next quickly grows as more time steps are 
taken). Generated particle tracks also allow for arbitrary 
initial conditions, a great advantage when exploring the 
behavior of coherent structures. The generated particle 
tracks solved the equation 

§ = u(r(i),t), (7) 

where r(t) is the position of the tracer particle at time t, 
using bicubic interpolation to approximate u at the par- 
ticle position and fourth-order Runga Kutta to perform 
the time integration. Typical time-spliced particle tra- 
jectories and generated trajectories are shown in Fig. |21 
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FIG. 3: Typical time-spliced real particle trajectory (□) and 
trajectory computed from dynamic velocity field data (•). 
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B. Coherent structure identification: the BPHK 
technique 

As stated in Section ^ we want to understand how the 
statistics of the Lagrangian dynamics depend on the exis- 
tence of coherent structures, such as vortices and saddles. 
We use here an Eulerian coherent structure census tech- 
niques to approximate the coherent structure regions in 
the flow using the interpolated velocity fields obtained 
through particle tracking. The BPHK criterion ^ ^ 
for the existence of a coherent structure is a constrained 
form of the Okubo- Weiss criterion [22l |2^ ■ For this rea- 
son, following the presentation given in |2(| , we start with 
a brief description of Okubo- Weiss. Assume that we are 
dealing with an inviscid 2D fluid. Given this assumption, 
the time evolution of the vorticity field is given by 



dt dt 1 dxi 



(8) 



The equation simply states that the vorticity following a 
fluid parcel is conserved in Euler fluids. Consider, for the 
same fluid, the evolution of the vorticity gradient, gi = 
d{U), also following the path of the fluid. The equation of 
motion for g is: 



dt 



(9) 



where the the velocity gradient tensor is given by Aij = 
diUj. 

Assume that the matrix A evolves slowly compared to 
the vorticity gradient g. Then one may consider A as es- 
sentially constant, and the vorticity gradient's evolution 
is determined locally by the velocity gradient. Since A is 
traceless by incompressibility, the local topology is deter- 
mined by |A|. In the case where |A| > 0, the eigenvalues 
of A are imaginary and the vorticity gradient is subject 
to a rotation. Otherwise, the eigenvalues are real and 
the vorticity gradient is subject to a strong strain, which 
enhances the vorticity gradient. Thus within the Okubo- 
Weiss approximation there are two types of topological 
structure within the flow, centers for which |A| S> and 
saddles for which |A| <C 0. 

The critical assumption in Okubo- Weiss, as pointed 
out by Basdevant and Philipovitch |20j, is that the ve- 
locity gradient matrix evolves slowly compared with the 
vorticity gradient. This assumption can be made more 
quantitative by taking the material derivative of Eq. 
This yields 



alA,, 



alt 2 



9j- 



dt 



3 dt 



(10) 



For the Okubo- Weiss assumption to be true, the second 
term must be much smaller than the third. That is 



< 1. 



(11) 




FIG. 4: BPHK coherent structure regions superimposed on 
the corresponding vorticity field. Ovals enclose centers (or 
vortices) and crosses denote the saddle positions. 



Without tracing the algebra (which is clearly enumerated 
in either [2(J or HJ), the above inequality can be re- 
written in terms of the eigenvalues X± of the Hessian of 
pressure, ify 



dfjp. Using these, Eq. ^] becomes 



R 



(A+ + A_)2 



< 1. 



(12) 



The implementation of BPHK is now fairly straight- 
forward. From the measured velocity fields evaluate the 
pressure field by inverting V 2 p = 2|A|. Obtain the 
pressure Hessian and evaluate the field R. Find points 
for which R < 1 corresponding to positions where the 
Okubo- Weiss assumption is valid. Now evaluate |A| at 
these positions and separate them into regions for which 
A | < (saddles) and |A| > (centers). These points 
comprise the coherent structures. The coherent structure 
identified using the BPHK criterion superimposed on the 
corresponding vorticity field is displayed in Fig. 



C. System characterization 

Under the driving conditions and for the magnetic con- 
figuration described in Sec. Ill Al the turbulence in the 
stratified layer system described above can be character- 
ized by several of its Eulerian ensemble-averaged statis- 
tics. The average energy per unit mass in the system 



E = u 2 . ms /2 was 8.4 enr/s (u rms = 4.1 cm/s) and 
the average enstrophy O = w 2 ms /2 where the vorticity 
uj = V x u was 51 s -2 (w rms = 10.1 s _1 ). The spectra of 
energy E(k) and enstrophy Q,(k) calculated using a Welch 
spatial window are shown in Fig. |SJ There is a build up 
of energy at large scales as expected for a system with 
an inverse cascade of energy. However, the spectra never 
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FIG. 5: The energy and enstrophy spectrum obtained from 
the electromagnetic layer for the driving conditions described 
in Sec.|n] Inset are the same spectra on a log-log plot. 



achieve the fc -5 / 3 for E(k) and fc 1 / 3 for fl(k) consistent 
with the existence of an energy inertial range. Also note 
that the spectra are significantly steeper in the enstrophy 
range than the predicted fc~ 3 , possibly due to the three 
dimensionality of the system at small scales. Using the 
spectra we can obtain the outer scale approximation of 
Lout ~ 6 cm and injection scale velocity, Ui n j — 1.8 cm/s 
and vorticity uji n j — 5.8 s _1 . 

The energy dissipated by viscosity and the frictional 
damping are given by e v = vuj^. ma and e Q = cra 2 ms 
respectively. Using the coefficients a and v given in 
Sec. lITAl yields e„ = 1.2 cm 2 /s 3 and e Q = 2.1 cm 2 /s 3 . 
Since the system is in an energetically steady state this 
implies that the energy injection rate e, n j = 3.3 cm 2 /s 3 . 
That more energy is lost to friction is consistent with an 
inverse cascade that moves energy predominantly away 
from small scales where viscous dissipation is efficient. 
One can also use the estimate e; ry « u^ n j/ri n j which 
yields a value of 2.9 cm 2 /s 3 , consistent with the above 
value. Similar results can be obtained for the enstro- 
phy injection and dissipation rates. For these we obtain 
f3 inj = 39.2 s" 3 , (3 a = 12.3 s~ 3 and /3„ = 26.9 s~ 3 . 

The bulk dissipation rates can be further character- 
ized by measuring the transfer of energy e(r) and en- 
strophy (3{r) through a length scale r using the filter 
approach reported in 01 • These functions are displayed 
in Fig. The flux directions are correct, predominantly 
downscale for enstrophy and upscale for energy. It should 
be stressed that the choice of filter used in the filter ap- 
proach, namely a Gaussian, has a tendency to smooth 
features; different filters would produce somewhat flatter 
regions (i.e. more inertial) in the energy and enstrophy 
transfer and sharper behavior around the energy injec- 
tion scale. Use of such filters, however, is prohibited by 
the existence of measurement boundaries. 

Finally, several length and time scales can be extracted 



from the reported coefficients. The Kolmogorov scale 
r] = (u 3 1 (3inj) 1/6 ~ 0.07 cm and the Taylor scale A = 
(E/Q) 1 / 2 » 0.41 cm, from which we can obtain the Tay- 
lor scale Reynolds number Re\ = \u rms jv = 112. Sim- 
ilarly the injection scale and outer scale Reynolds num- 
bers are Re = ri n jUi n j/v — 303 and Re = L u rms /v = 
2150. We should note that -q is small, much less than 
the salt layer depth of 0.3 cm. We therefore expect 
3D effects to become important before the Kolmogorov 
scale is reached. Characteristic times that one can define 
for the turbulent state are the large-eddy turnover time 
t = L /u rms = 1.5 s, the injection scale eddy turnover 
time Ti n j — ri n j/ Ui n j = 1 s and a characteristic shear 
time t = n~ 1/2 = 0.14 s. 



III. LAGRANGIAN STATISTICS 

Before presenting results involving coherent structures, 
Lagrangian single- and two-particle statistical quantities 
for the full data fields are considered. In this section 
both real and generated particle trajectories are used in 
the analysis. 



A. Single trajectory statistics: Self-diffusion 

Single trajectory statistics depend on how the time in- 
crement (t— to) compares with the Lagrangian correlation 
timeT L =/~(u(r(i ),i )-u(r( S ), S ))d S /(|u| 2 }. Figure[7| 
shows the velocity autocorrelation function for both La- 
grangian trajectories and Eulerian points. The latter is 
displayed to contrast the 2D turbulence of the stratified 
layer with 3D turbulence results where te < tl [te has 
the same definition as tl except that the "trajectory" is 
a constant point in space). Performing the integration 




FIG. 6: The scale-to-scale transfer of energy and enstrophy 
obtained from the electromagnetic layer for the driving con- 
ditions described in Sec. EI 
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FIG. 7: The Lagrangian and Eulerian velocity correlation 
function obtained using real particle tracks from the stratified 
layer. Using these functions the Lagrangian correlation time 
(see SecQfor definition ) is found to be tl = 0.26 s. Similarly 
the Eulerian correlation time is te = 0.95 s. 



and normalization yields t l = 0.26 s and te = 0.95 s. 
The ratio te/tl — 3.7 is consistent with earlier numer- 
ical simulations of 2D turbulence which found the ratio 
to have values between 3 and 4 The oscillations in 
the Eulerian autocorrelation function, corresponding to 
the 0.5 Hz square wave in the electromagnetic driving 
current, are not apparent in the Lagrangian autocorrela- 
tion function because the particle trajectories randomly 
sample spatial points, averaging out the oscillations. 

The second moment of single particle displacements 
(|5r| 2 ) is displayed in Fig. [5] as a function of time dis- 
placement (t — to) for both real and generated trajecto- 
ries. The agreement between the two is fair, the largest 
difference being at large scales where imposing artificial 
boundaries on the generated particles limits the max- 
imum measurable displacement and reduces the aver- 
age. There is a slight disagreement at approximately 
2t l where the generated data tends to overshoot the real 
data. This may be statistical bias in the particle tracking 
algorithm (splicing tracks over smooth regions is easier 
than over heavily strained regions) or could be an indi- 
cation of inertial behavior of the real particles. 

The trend in (|<5r| 2 ) for it—to) < Lis power-law scaling 
with an exponent ~ 2 as expected in the ballistic regime 
of motion. The amplitude is 14 cm 2 /s 2 , fairly close to the 
expected value of w 2 ms = 17 cm 2 /s 2 . For (t — 1 ) > L the 
exact exponent is sensitive to the boundary, but a power 
law fit to the data without imposed boundaries yields an 
exponent of about 1.3 in agreement with earlier numer- 
ical results 0. Finally, no Brownian diffusive regime is 
reached in the current data set. To see this regime, data 
from a significantly larger spatial area would be required. 



B. Two trajectory statistics: The Richardson 
regime 

As discussed in Sec. the expected behavior of the two 
trajectory Richardson statistics in an inverse cascade in- 
ertial range is (|R| 2 ) ~ St 3 . In Fig. [3]the two-trajectory 
statistics for (|R| 2 ) are displayed for several initial dis- 
placements |Ro|. A very similar set of curves were ob- 
tained from numerical simulations of 2D turbulence p|. 
Inset are the same plots compensated by (t — to) 3, . It 
should also be noted that the expected exponential small 
time behavior (not shown) is not a clean exponential for 
the generated data. Moreover, the dependence of (|R| 2 ) 
is markedly different for real and generated data at small 
time, though long time behavior is identical (real data 
tends to show more of an t 2 scaling at small time possi- 
bly due to the real particles inertia). Due to statistical 
constraints, the generated data will be used in subsequent 
analysis. 

For an initial displacement of 0.075 cm, there seems to 
be a nice 8t 3 range from 0.3 s to 1.0 s. That the scaling 
depends on the initial separation and that the spatial 
scale where scaling is apparently observed is below the 
injection scale are two reasons to be cautious about inter- 
preting this scaling as indicative of the predicted Richard- 
son result 0. Further, one expects crossover effects 
to play an important role in the precise determination 
of scaling for moderate Reynolds number flows. Details 
of Richardson scaling and comparisons with results from 
other experiments|l3|| will be presented elsewhere. Here 
the effects of conditioning the separation on whether the 
initial positions of the particles is within saddles or cen- 
ters is considered. Unless otherwise noted the initial dis- 
tance of |Rq| = 0.75 cm will be used for the rest of the 



10 2 P j 




5t/x L 



FIG. 8: Single-particle Lagrangian displacement (|<5r| 2 ) for 
real and generated trajectories. With the generated trajecto- 
ries several different boundary sizes were utilized to illustrate 
the effect of finite data extent on the results. 
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FIG. 9: The mean square relative displacement (|R| 2 ) of 
two-particle trajectories in the flow for several different initial 
separations i?o- All of the plots were created using generated 
particle trajectories. Inset are the same plots normalized by 
the Richardson predicted behavior of St 3 . 
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FIG. 10: The conditional Lagrangian velocity correlation 
function obtained using generated particle tracks from the 
stratified layer. Dashed and dash-dot lines correspond to con- 
ditioning upon the type of structure the trajectory is within 
at time to. The solid line is the unconditioned average. In- 
set is the probability distribution for the topological quantity 



analysis. 



IV. EFFECTS OF COHERENT STRUCTURES 

In this section the statistical quantities measured in 
Sec. IIIII for ensembles of particles starting at any given 
point in the flow are conditioned on whether the initial 
particle position r(irj) is within a saddle or center as iden- 
tified using the BPHK technique described in Sec. [n] All 
results presented in this section were obtained from gen- 
erated data rather than real particle tracks so that statis- 
tical convergence could be achieved. In general one would 
like to separate the statistics on the basis of how long that 
particle was within the coherent structure and how long 
that coherent structure existed. This separation would 
be particularly important for the case of saddles which do 
not have a closed group of core particles associated with 
them. Unfortunately, such conditions further constrain 
statistics that have already been pushed to their limit, 
and thus, only the approach based on initial conditions 
is possible. 



A. Conditional single trajectory statistics 

As in Sec. IIII Al the discussion of single trajectory 
statistics starts with the single point velocity auto- 
correlation. The correlations are displayed in Fig. ^3 
where the average () has been conditioned on the type of 
coherent structure the particle is within at time to', the 
average over the full data field is shown for comparison. 

For both saddles and centers, the magnitude of the cor- 



relation function at St = is suppressed compared with 
the background flow. The suppression is caused by the 
self-organization of coherent structures which are rela- 
tively quiescent compared with the background flow, and 
thus they have lower RMS fluctuations. One might be 
surprised that this effect is stronger for saddles than for 
centers. The difference is due to a bias in the population 
statistics of saddles and centers |3l l25j : it is more prob- 
able to find a center of significant strength in the flow 
than a strong saddle. This argument is supported by 
considering the probability distribution, P(|A|) for the 
topological quantity |A| defined in Sec. Ill Bl and shown 
in the inset of Fig. ^5] There is a larger probability for 
|A| 3> than |A| <C 0. For more discussion about why 
this must be the case see |25| . 

For both saddles and centers, when 8t 1 s this initial 
difference has, for the most part, disappeared and the 
correlation functions collapse to the correlation function 
of the full statistics. For centers, the correlation func- 
tion oscillates with a period of about 0.6 s, fairly close 
to the injection scale eddy rotation time, Tj n j = 1.0 s, 
which is expected owing to the periodic rotation of the 
centers. The steep fall off in the amplitude of the center- 
correlation occurs for a number of reasons: the stripping 
of particles from the exterior of vortices, the random walk 
of vortices through the flow, and a finite vortex lifetime. 
In spite of these effects, weak oscillations in the corre- 
lation function continue out as far as 2.4 s, indicating 
some long lived structures in the flow. Using the center- 
correlation function, the eddy turnover time is defined to 
be r e = 0.6 s, which is close to the dimensionally esti- 
mated value of the injection-scale eddy turnover time of 
1.0 s. 

The behavior of the correlation function for saddles is 
somewhat surprising in that it initially increases. This 
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increase comes about because particles starting within 
a saddle are not confined to the saddle, as they are in 
the cases of centers, but are leaving the saddle to the 
background. As already indicated, the background has 
stronger fluctuations when compared with the relatively 
calm region within the saddle, thus the increase. The 
crossing point of the saddle-correlation with the full cor- 
relation at about 0.25 s is interpreted to be the average 
time a particle remains within a saddle structure. After 
this time, out to about r e , the saddle correlation exceeds 
that of the full statistics. This crossover indicates that 
the particles, once having left the saddle, continue to 
remember their interaction with the structure for quite 
some time, effectively surfing the saddle's unstable man- 
ifold. 

Given the observations obtained from the Lagrangian 
velocity correlation, the conditional single trajectory 
statistics, displayed along with the unconditioned statis- 
tics in Fig. 1111 can be interpreted. One might expect 
that at small times, corresponding to being within sad- 
dle structures, the self diffusion would be enhanced but 
at very small times self diffusion within saddles is smaller 
than for any other point in the flow. This slow behav- 
ior can be explained by the suppressed RMS fluctuations 
within saddles resulting from a bias in the population 
statistics of |A|. A similar argument explains the sup- 
pression of self diffusion within centers at early times. 
The expected enhancement in self-diffusion within sad- 
dles does not occur until times greater than (t—to) = 0.25 
s, where the earlier conditional Lagrangian velocity cor- 
relation crosses the full Lagrangian correlation function. 
The crossing time is where most of the trajectories have 
left the saddle. The conclusion, then, is that saddles en- 
hance diffusion, but non-locally. Only when the particle 
escapes the confines of the saddle structure and begins 
to surf the unstable manifold is the enhancement felt. 
At long times, about 2 s, the trajectory begins to forget 
both the saddle and unstable manifold and the condi- 
tional statistics approach the full statistics. 

The self diffusion for particles beginning within cen- 
ters behaves more intuitively. The suppression at small 
times was explained by bias in the statistics of |A|. At 
intermediate times, smaller than the eddy rotation time, 
self-diffusion continues to be suppressed because the tra- 
jectories cannot cross the center manifold and continue 
to circulate about the vortex. Finally, at long times, as 
the centers dissolve, the conditional statistics approach 
the full statistics. Overall, the effect of different initial 
conditions is not large, accounting for at most a factor of 
2 between centers and saddles. 



B. Conditional two-trajectory statistics 

The results of conditioning the average (|R| 2 ) upon the 
condition that a trajectory starts within a center or sad- 
dle at to are displayed in Fig.^|. In these two-trajectory 
statistics the expected enhanced diffusion due to saddles 
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FIG. 11: Single-particle Lagrangian displacement (|<5r| 2 ) for 
generated trajectories conditioned upon the type of structure 
the trajectory is within at time to. The solid line represents 
the unconditioned average. 



is realized. From the earliest times to times as long as 2 
s there is an enhancement in the relative diffusion with 
a peak difference at about 0.25 s. From Sec. IIV Al rc- 
call that this earlier time is roughly the time it takes 
for the particles to leave the saddle structure. Thus, for 
relative diffusion, in stark contrast to the self-diffusion 
presented before, saddles behave locally. Once the max- 
imum amount of enhancement occurs, the saddle statis- 
tics decay back to the full statistics as the particle pairs 
forget their origin within the saddle. The decay is slow 
most likely because some particle pairs initially straddle 
the stable manifold and thus exit the saddle on opposite 
sides and surf opposing unstable manifolds. At around 2 
s, however, the statistics collapse back to the full statis- 
tics, much as in the case of self-diffusion, because the 
trajectories begin to forget their origins. 

The relative diffusion of particles within centers also 
behaves, more or less, as expected. The diffusion is sup- 
pressed at small times and slowly approaches the full 
statistics as the vortices are destroyed and the trapped 
trajectories are released to evolve in the background. One 
might be surprised that the trapping of vortices does not 
suppress the diffusion with respect to the full statistics 
more than the slight amount displayed here. To resolve 
this issue, consider the types of structures, other than 
saddles, that exist within the flow. Whatever their form, 
for example jets, one expects that by definition |A| is 
small, i.e., there are small velocity gradients. Thus, in 
these regions, a pair of particles would evolve together, 
only separating quickly when they encounter a saddle. 

To better visualize the implications of the above re- 
sults on relative diffusion 10 5 trajectories were generated 
from initial positions within a center region and within 
a saddle. The results of tracking these trajectories over 



9 




(t-g/T. . 

v (y inj 



competing scalings at low to moderate Reynolds number 
and of boundaries. The intermediate time behavior de- 
pends on the initial starting separation Ro with decreas- 
ing slope for large Ro. The existence of a Richardson 
scaling range with (R 2 ) oc (t — to) 3 is unresolved and will 
be presented in detail elsewhere. At short times the sep- 
aration is sensitive to starting conditions, i.e., whether 
the initial mean position was within a saddle or cen- 
ter. As expected, separation increased more rapidly at 
small times when the particles started in a saddle and 
less rapidly when they were within a center. 

The most interesting observation for the separated 
statistics is the time scales at which the structures, in 
particular the saddles, have their greatest effect. For the 
two-point statistics the greatest effect of structures is felt 
at time scales for which the particle is within the struc- 



FIG. 12: Two-particle Lagrangian displacement (|R| 2 ) for 
generated trajectories conditioned upon the type of structure 
the trajectory is within at time to- The solid line represents 
the unconditioned average. 

time are displayed in Fig. 1131 One can see from the 
figure that the initial spread for trajectories beginning 
within the saddles is far greater than for the correspond- 
ing center region. It should be stressed that only a few of 
the starting particles within the saddle account for this 
large spread, the majority continue to exist in the dense 
grouping of particles in the left of the image. This dense 
grouping weights statistics, such as those shown in Fig. 
1131 so that they only echo a fraction of the magnitude of 
the spread. 



V. CONCLUSIONS 

The Lagrangian approach of particle diffusion was ap- 
plied to the quasi-2D turbulence of a thin stratified 
layer. Simultaneous cascades of enstrophy and energy 
develop although neither is strictly inertial despite larger 
Reynolds number and smaller damping than in other 
stratified layer experiments ^(|. Single particle diffu- 
sion follows the predictions of Taylor at small times that 
(5r 2 ) oc (|u| 2 )<5t 2 and is consistent with the scaling of 
(t — to) found from numerical simulations for times 
larger than the Lagrangian velocity correlation time[j|. 
The expected long-time behavior of diffusive scaling like 
t — to is not observed owing to the finite size of the sys- 
tem. By separating the statistics based on the initial 
position within a topological saddle or center, one can 
understand the influence of these structures on particle 
diffusion. The overall effect is small with the maximum 
difference near Ti n j where center dispersion is about half 
that of particles in the rest of the field. 

The behavior of two-particle separations is more com- 
plex, reflecting the competing effects of crossover between 
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FIG. 13: The time evolution of an "ink drop" of generated 
particle tracks. The left column represents the behavior of an 
ink drop that initially seeds a center region. The right column 
is the same for a saddle region. 
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ture. The single point statistics, on the other hand, are 
non-local, that is the structures greatest effect on these 
is at times when the particle has exited the structure. 
This difference might be explained by the existence of an 
unaccounted for structure, e.g. a jet, which lurks on the 
outside of saddles. When a particle is advected out of 
the saddle, it is swept up in a jet and advected rapidly 
away from its starting point, which explains the enhanced 
self-diffusion at intermediate times. If two particles are 
advected into a jet they are carried along with little rel- 
ative displacement, thus the two-point statistics are not 
enhanced. For this explanation to be compelling, the 
notion of a jet needs to be quantitatively defined. 

The potential for the Lagrangian approach to the anal- 
ysis of turbulence has barely been scratched by these 
preliminary studies. Multiple-point particle tracking can 
be implemented easily and would provide a better mea- 
sure of turbulence transfer processes. Another approach 
would be to measure interscale transfer flux in the mov- 
ing frame of a group of particles that form a particular 



topological structure, e.g., saddles or centers. One can 
compute the stretching fields of the flow to probe the mix- 
ing properties of the turbulent system |2a|. There are 
many more analysis ideas that follow along these lines 
and that take advantage of being able to measure the 
dynamics of fluid quantities moving with a fluid parcel. 
The opportunities are numerous and will help elucidate 
the mechanisms for turbulence and transport in 2D tur- 
bulence. Further, the concepts and analysis that can be 
applied to 2D data will set the stage for similar applica- 
tions to 3D data when it becomes available. 
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